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We consider two-dimensional dipolar bosonic gas with dipoles oriented perpendicularly to the 
plane in a weak random potential. We investigate analytically and numerically the condensate de¬ 
pletion, the one-body density-matrix, the ground state energy, the sound velocity and the superfiuid 
fraction. Concentrating on the regime where a rotonlike excitation spectrum forms, our results show 
that the superfluidity disappears below a critical value of disorder strength yielding the transition 
to a non-trivial quantum regime. 
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Recently, ultracold dipolar gases in two-dimentional 
(2D) geometry have been the subject of intense exper¬ 
imental and theoretical investigations [H-Q- What ren¬ 
ders such systems particularly intriguing is the presence 
of the low-lying roton minimum in the excitation spec¬ 
trum i0’ and the possibility of the crystallization of 
solid bubble into a lattice superstructure, resulting in 
a global supersolid phase i i 0. However, such su¬ 
persolids require a dense regime with several particles 
within the interaction range, which can be difficult to 
achieve. The same holds for supersolids discussed for 2D 
dipolar Bose gases [l^ near the gas-solid phase transi¬ 
tion [ij, [T^ . It was found also that this state appears in 
B-t-J model of two-component bosons as a result of the 
long-range DDI (l3j |. 

A more complicated situation arises if these phenom¬ 
ena are studied in a random environment. Bose gases 
in the random medium attracted a great deal of interest 
because it connects two central ideas of the condensed 
matter theory: the Bose-Einstein condensation (BEC) 
and localization. This latter occurs both for fermions and 
bosons 0 and results in the existence of an insulating 
phase called Bose glass 0. The existence of a phase 
characterized by simultaneous glassiness and superfluid¬ 
ity and thus constitutes a glassy counterpart to the su¬ 
persolid phase 0-01, was first observed in numerical 
Quantum Monte Carlo simulations of solid "^He samples 
by Boninsegni etaZ0|. This phase which is featured (at 
the same time) by superfluidity and a metastable amor¬ 
phous structure called ” superglass”. Superglass may also 
have realizations in interactin g b osons at very low tem¬ 
peratures and high density [1,0. In this context, Zam- 
poni etal have shown that quantum fluctuations can 
stabilize the superglass phase in a self-disordered envi¬ 
ronment induced by geometrical frustraion. 

On the other hand, many interesting works on BEC 
and superfluidity (See, e.g., [2ll - [2l| l have been reported 
for disordered cold atomic gases with pure contact inter¬ 
action in continuum, the study on dipolar boson systems 
is still inadequate. In the present paper, for the first 
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time to our knowledge, we investigate the properties of a 
quasi-2D homogeneous dipolar Bose gas in a weak ran¬ 
dom potential with delta correlated disorder. Our study 
is based on the Bogoliubov approach, this method which 
marked an important step towards a quantitative descrip¬ 
tion of dirty dipolar Bose systems 0, 0|, allows for 
accurate determination of the condensed depletion, one- 
body density matrix, ground state energy, sound velocity 
and superfluid fraction. It is found that the presence of a 
disordered potential in the regime where the roton devel¬ 
ops in the excitations spectrum strongly enhances fluctu¬ 
ations and thermodynamics quantities. We demonstrate 
that both BEC and superfluidity are depressed due to 
the competition between disorder and DDI yielding, the 
transition to an unusual quantum phase. 

We consider a dilute Bose-condensed gas of dipolar 
bosons in an external random potentials. These parti¬ 
cles can be confined to quasi-2D, by means of an external 
harmonic potential in the direction perpendicular to the 
motion (pancake geometry) and all dipoles are aligned 
perpendicularly to the plane of their translational mo¬ 
tion, by means of a strong electric (or magnetic) field. In 
this quasi-2D geometry, at large interparticle separations 
r the interaction potential is V(r) = dP/r^ = fpTif/mr^, 
with d being the dipole moment, m the particle mass, 
and r* = md?/fP the characteristic dipole-dipole dis¬ 
tance. The disorder potential is described by vanishing 
ensemble averages {U (r)) = 0 and a finite correlation of 
the form {U(r)U{r')) = R{r,r'). 

In the ultracold limit where the particle momenta sat¬ 
isfy the inequality fcr, <C 1, the scattering amplitude is 
given by (see e.g. i) 

f{lk')=g{l-C\k-k'\), (1) 


where the 2D short-range coupling constant is g = 
gsD / V^lo with /q = u) is the confinement fre¬ 

quency and C = 2TTtPr^,/mg = 2TTd^/g. EmplOTing this 
result in the secondly quantized Hamiltonian (6|, we ob¬ 
tain 
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H-^^Ekdtd^+ ( 2 ) 

k k,p 

k,q,p 


The noncondensed density is defined as h = J2k(^\^k)■ 
Then invoking for the operators aj: the preceding Bo- 
goliubov transformation, setting = dj:pN}^ and 

putting the rest of the expectation values equal to zero, 
where Nk = [exp(efe/r) — 1]“^ are occupation numbers 
for the excitations. Using the fact that 2N{x) + 1 = 
coth(a;/2) [2^, we obtain: 


where S is the surface area, Ek = Ti^k'^/2m, and at, d^ 
are the creation and annihilation operators of particles. 
At zero temperature there is a true BEG in 2D, and we 
may use the standard Bogoliubov approach. Assuming 
the weakly interacting regime where mg/2nh'^ 1 and 
r* <C with ^ = h/firing being the healing length. We 
may reduce the Hamiltonian Q to a bilinear form, using 
the Bogoliubov transformation 21| a^: = Ukb^ — vp - — 


/3g, where St and S^ are operators of elementary exci¬ 
tations. The Bogoliubov function s Uk, V k are expres sed 
in a standard way: Uk,Vk = {\/ek/Ek ± ^/Epek)l2, 
= pnJSUkEk/e^, and the Bogoluibov excitation en¬ 
ergy is given by et = \/El + 2ngEk{l - Ck). 

To zero order the chemical potential is g = ng. If 
C/e < V8/3, Ek is a monotonic function of k. However, 
it shows a roton-maxon structure for the constant C in 
the interval a/S/S < C/£< 1. It is then convenient to 
represent Sk in the form [^: 


£k = 


fPk 

2m 


\J{k- krY + k\, 


(3) 


where kr = 2ClPp and fcA = \/— fc/. If the ro- 
ton is close to zero, then kr is the position of the roton, 

and A = ?Pkrk^/2m = 2ngC\Jmng/tP — C'^{mng/Tp)'^, 
is the height of the roton minimum. At C/^ = 1, the ro¬ 
ton minimum touches zero and for C/^ > 1, the uniform 
Bose condensate becomes dynamically unstable and the 
uniform superfluid is no longer the ground state. 
Importantaly, the spectrum (|3]) is independent of the ran¬ 
dom potential. This independence holds in fact only in 
zeroth order in perturbation theory; conversely, higher 
order calculations render the spectrum dependent on the 
random potential due to the contribution of the anoma¬ 
lous terms {d_-^d-p). 

The diagonal form of the Hamiltonian of the dirty 
dipolar Bose gas can be written in the usual form 
H = E + where E = Eq + SE + Er with 

Eq = SgrP/2 and SE = i '^[sk — Ek — ng{l — Ck)] be- 

k 

ing the ground-state energy correction due to qunatum 
fluctuations. 


Er 


r ^k 


y^nRk 

k 



(4) 


gives the correction to the ground-state energy due to the 
external random potential. 


■'= 2 / 1 : 


Ek -I- ng{l - Ck) 


£k 


coth 


\2TJ 


- 1 


where 
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riR, 

(5) 

( 6 ) 


is the condensate depletion due to the external random 
potential. 

In order to investigate in a simple way how the random 
potential affects the behavior of the system, we will often 
make the white-noise assumption in which the external 
potential is described by a single parameter i?(r, U) = 
RoS{f,r') [ 2 ^, where Rq denotes the disorder strength 
which has dimension (energy) ^ x (length)^. 

Let us now asume that the roton is close to zero and the 
roton energy is A ^ nt;, we have the coSicient C close to 
and kr ~ 2/^. Then, using Eq.®, for the contribution 
of momenta near the roton minimum at T = 0, we obtain: 


^ /‘^ngP 

n \ A ) 


(7) 


where R = Rq jng^ is a dimensionless disorder strength. 
The noncondensed density is calculated via m as 


n 


mg 

ttTtP 




A < ng, (8) 


where = ^JPPPJPmg. 

The leading term in Eq.® which comes from the quan¬ 
tum fluctuations was first obtained in our recent work , 
while the second term which grows faster than the first 
one, represents the disorder correction to the conden¬ 
sate depletion. The condensed fraction can be written as 
ric/n = 1 — n' jn. 

To check the result of Eq.® , we solve Eqs.® and ® 
numerically using Monte Carlo method. FigurelT] shows 
that in the absence of the random external potential i.e. 
i? = 0, the noncondensed fraction grows logarithmically 
(see Eq. ®) when the the roton energy A goe s to zero 
yielding the transition to a supersolid state @,ii. It 
is worth stressing that in the context of the liquid he¬ 
lium, the position of the roton minimum influences also 
the phenomenon of superfluidity 27| . 

In the presence of the external disordered potential and 
for C ~ the condensate depletion becomes more signif¬ 
icant and diverges when A = 0 indicating the transition 
to a novel quantum phase in dilute 2D dipolar bosons. 
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FIG. 1. Quantum depletion of the condensate from Eq. 
as a function of G/^ for — 0.01. Clean dipolar 

BECidashed line {R = 0). Dirty dipolar BEG: solid line {R = 
0.5). 


FIG. 2. One-body density correlation function for rngjA-Kh'^ = 
0.01. Clean dipolar BECidashed line {R — 0). Dirty dipolar 
BEG: solid line [R = 0.5). 


The same factors of Eq.® appear in the one-body den¬ 
sity matrix gi(f) = (Tf(r)'F(O)), where ^(r) is the field 
operator. The correction due to the disorder effects to 
the correlation function is g^{f) = f 
Assuming that the roton minimum reaches zero and tak¬ 
ing into account only the contribution of momenta near 
this minimum we have for A <C 


Equation (ED shows that grows linearly with ng/A, 
and has a negative value which leads to reduce the total 
energy of the system. 

The correction to the ground-state energy due to quan¬ 
tum fluctuations can be given as 


Eo 


2m g 
ttTi^ 



A <C ng. (12) 


9i{r) ^ ^ rr^ 
n 



-R 



M2r/0, 


(9) 

where Jo{x) is the Bessel function. 

The numerical calculation of the one-body density matrix 
well agrees with the analytical result obtained in Eq Q 
and shows that gi (r) decays at large distance when C ~ ^ 
as is depicted in Fig® This signals the non-existence 
of the off-diagonal long-range order (i.e., BEC) in the 
disordered 2D dipolar bosons. 

As we see from Eq. ® , for the roton minimum close to 
zero a small condensate depletion requires the inequality 



( 10 ) 


We thus conclude that at T = 0, the validity of the Bo- 
goliubov approach is guaranteed by the presence of the 
small parameter ED- 

However, the situation changes in the calculation of 
the correction to the ground-state energy due to the ex¬ 
ternal random potential. When the roton minimum is 
approaching to zero, we get from (|T]) 

Er 2mg (2ng\ 


The disorder correction to the chemical 
calculated esealy through dEn/dN 


9‘R 




V A / 


potential can be 

A <C ng. (13) 


Note that quantum fluctuations corrections to the chem¬ 
ical potential had already obtained in our recent paper 


One can also show that the shift of the sound velocity is 
consistent with the change in the compressibility mc^ = 
ndg/dn and is given by 


4 


^sO 


TTh^ 





(14) 


where c^o = ygjm is the zeroth order sound velocity. 
The second and the third terms originate from quantum 
fluctuations while the last term comes from the disorder 
contribution. 

On the other hand, in an infinite uniform 2D fluid ther¬ 
mal fluctuations at any nonzero temperature are strong 
enough to destroy the fully ordered state associated with 
BEC, but are not strong enough to suppress superfluidity 
in an interacting system at low, but non-zero tempera¬ 
tures. The presence of residual ” quasi-long-range” order 
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at low temperatures leads to an interesting interplay be¬ 
tween superfluidity and condensation in all experimen¬ 
tally relevant finite-size systems. In this quasiconden¬ 
sate, the phase coherence governs only regimes of a size 
smaller than the size of the condensate, characterized by 
the coherence length 4[2i,[33. Thermodynamic proper¬ 
ties, excitations, and correlation properties on a distance 
scale smaller than 4 are the same as in the case of a true 
BEC. Upon utilizing the previous definitions we find that 
the correction to the condensate depletion, the correla¬ 
tion function and thermodynamic quantities due to ther¬ 
mal flucatuations is given by the factor (2mg/?i^)T/A 

a. 

The superfluid fraction n^/n can be found from the 
normal fraction rin/n which is determined by the trans¬ 
verse current-current correlator nsjn = 1 — nnjn. We 
apply a Galilean boost with the total momentum of the 
moving system Pv = P -I- mvN, where P = kata^: 
and V is the liquid velocity. In d-dimensional case, the 
superfluid fraction reads 


!4. - 1 _ ^ 

n dTn 


d‘^k 

(27r)‘^ 


Ek 


4sinh4efe/2r) 


-f 


nRoE? 


coth 


\2TJ 


(15) 


At very low temperature we can put coth(efe/2T) = 
2Tjek- Thus, Eq. (IT5|) reduces to 


!4. - 1 _ 

n d n 


2 

f d’^k 

Ek 

dTnJ 

[2'k)’^ 

.4sinh4efe/2T). 


• (16) 


Interestingly, the ratio between the normal fluid density 
and the corresponding condensate depletion increases to 
2 in 2D and to 4 in ID, in contrast to the familiar 4/3 
in 3D geometry obtained earlier in [2l|, l^. Another 
important remark is that the superfluid fraction m is 
no longer a tensorial quan tity as in the case of a 3D dirty 
dipolar Bose gas since the dipoles are assumed 

to be perpendicular to the plane. However, if the dipoles 
would be tilted slightly, superfluidity would acquire an 
anisotropy and thus becomes a tensorial quantity. 

Assuming now that the roton minimum is close to zero 
and A ^ T, then the momenta near the roton minimum 
are the most important, and the use of Eq. © yields: 


n- 


mg (2ng\ 2mg T 

- 


A' 


(17) 


At r = 0 and in the absence of disorder, the super¬ 
fluid fraction is equal to unity. Another interesting con¬ 
sequence of the above results is that for any value of 
(ng/A) there exists a critical strength of disorder Rc = 
(4/7r) ln(2ng/'/A)/ (2ng/A)^ for which rig/n < ric/n. For 
dysprosium atoms at 2D densities ^ 10® cm“®, one has 
r* ~ 200 A, ng ^5 nK, mg/2'iT‘li? ~ 0.025 and A should 
be above 2 nK. Therefore, the precedent criterion is sat¬ 
isfied for i?c = 8 X 10“'*, which means that ng can be 
smaller than ric only for a very weak disorder i.e. i?c <C 1. 
One can conclude that the Bogoliubov approach should 
satisfy the conditions R < Rc- However, it is not clear 
whether these results still apply for R > Rc in a range of 
densities where the difference between rig/n and Uc/n can 
be significant. The response to these questions requires 
either a non-perturbative scheme or numerical Quantum 
Monte Carlo simulations, which are out of the scope of 
the present work. 


In conclusion, the impact of a weak random poten¬ 
tial on BEC and on superfluidity in a dilute quasi-2D 
dipolar Bose gas is studied with a combined numerical 
and analytical approach. Our analysis signifies a more 
pronounced effect of disorder in such a system when the 
roton is approching zero with enhancing quantum fluctu¬ 
ations, one-body density-matrix, equation of state, sound 
velocity and depleting superfluid density. Furthermore, 
we have reproduced the expression of the condensate 
fluctuations and thermodynamic quantities obtained in 
our recent work @ in the absence of the disorder po¬ 
tential. Our results represent a starting point for the 
analysis of collective modes of homogeneous or trapped 
2D dirty dipolar BECs. We have found, in addition, the 
criterion of applicability of the Bogoliubov approach for 
these systems. An important step for a future work is to 
show whether these disordered quasi-2D dipolar systems 
promote a stable superglass phase, where superflow and 
glassy density localization coexist in equilibrium without 
exhibiting phase separation. A qualitative study of such 
a phase diagram necessitates a more sofisticated scheme 
with a non-perturbative solution. 


We are grateful to Nikolay Prokofiev and Axel Pel- 
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manuscript. 
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